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Abstract: We study the cosmological perturbations in teleparallel dark energy models 
in which there is a dynamical scalar field with a non-minimal coupling to gravity. We find 
that the propagating degrees of freedom are the same as in quintessence cosmology despite 
that variables of the perturbed vierbein field are greater than those in metric theories. The 
resulting growth evolution shows that gravitational interactions are enhanced during the 
unique tracker evolution of teleparallel dark energy models. 
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1 Introduction 

The dark energy mystery has received much attention in the context beyond simply embed- 
ding the cosmological constant in general relativity (GR). To avoid a gigantic fine-tuning 
issue when the cosmological constant is attributed to the vacuum energy density, a canon- 
ical scalar field is introduced with the dynamical behavior known as quintessence [1-6]. 
Further modification of quintessence has been considered by allowing a non-minimal cou- 
pling between the scalar field and gravity [7-15] as well as either using the "wrong" sign 
phantom field [16-21], or combining a non-canonical extension such as the scalar-tensor 
theories [22, 23]. 

Recently, the "teleparallel" description of GR [24, 25] reveals an interesting aspect 
of gravity to explain the dark energy origin. The late-time cosmic acceleration has been 
seen from the generalized action of teleparallel equivalence of GR (TEGR) [26, 27], named 
as f{T) gravity where T denotes the gravity Lagrangian of TEGR. Although f(T) grav- 
ity shows a clear reference to f(R) gravity, it provides no higher than second-order field 
equations, and no trivial Einstein frame avalible through conformal transformations [28]. 
On the other hand, teleparallel gravity incorporating with a dynamical scalar field de- 
scribes a different scenario to drive the accelerated expansion of the Universe. This class of 
teleparallel dark energy models [29-31] lies in the framework of TEGR where the simplest 
quintessence model [1] is recovered if the scalar field is minimally coupled to gravity. Given 
that the torsion scalar T differs from the Lagrangian density of GR by a total derivative, 
the presence of a non-minimal coupling between the scalar field and gravity leads to in- 
teresting cosmological behaviors [1, 7, 16], while teleparallel cosmology shows compelling 
results with current observations [32]. 
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However, it is noteworthy that modified teleparallel gravity theories, including both 
f(T) gravity and teleparallel dark energy models, do not respect to the Lorentz invariance in 
the local tangent frame. The Lorentz violation in the teleparallel formalism beyond TEGR 
inevitably introduces extra degrees of freedom, which are unfamiliar in the framework of 
GR. This fact is manifestly provided by the field equations, which address the dynamics 
of all 16 components of the vierbein field [33, 34]. Under the homogeneous and isotropic 
principles of the background cosmology, those extra degrees of freedom merely involve at 
perturbation level, but have significant effects on the structure formulation in the Universe 
[35]. 

In the present work, we study the cosmological perturbations in teleparallel dark energy 
models via the perturbed vierbein field, which contains some variables with no reference 
in metric perturbations [36]. Nevertheless, in the scalar perturbations, we find that only 
one new degree of freedom involves at linear level, and we obtain one implicit constraint 
equation to the scalar mode variables from the asymmetric part of the field equations. As a 
result, the new scalar mode is resolved with a mere algebraic relation to other scalar modes, 
and thus the propagating degrees of freedom in teleparallel cosmology do not increase. 

We proceed the density perturbations in sub-horizon scales to study the growth be- 
haviors of some specific teleparallel dark energy models. For the simplest tracker solution 
corresponding to the non-minimal gravity-field coupling [37] (hereafter we denote as "NGF 
tracker"), we compare the growth rate from the effective gravitational constant to the 
GR case under an identical evolution background. We find that the NGF tracker always 
dominates, even in the presence of potentials of the scalar field, if the coupling between 
gravity and the scalar field is large to the order of unity. In the broader case with a smaller 
non-minimal coupling, the cosmic acceleration can be driven by the potential or the NGF 
tracker, depending on the initial composite differences of dark energy. The resulting density 
perturbation growth indicates that gravitational interactions are stronger if the evolution 
follows the NGF tracker rather than the usual potential-driven cosmic acceleration case. 

This paper is organized as follows. In Sec. II, we review teleparallel dark energy 
theories. In Sec. Ill, we demonstrate the perturbation equations and apply the results to 
study the matter growth in Sec. IV. Finally, conclusions are given in Sec. V. 



2 Teleparallel Dark Energy 

The teleparallel dark energy model is given by [29]: 

S-- 



d 4 xe 



(2.1) 



where e = det(e^) = y/—g with k 2 = 8irG, £ is the non-minimal coupling parameter, and 



rp jllT 1 [IVr-pp ^rpflU rpp rpp rpV p 



(2.2) 



is the torsion scalar of TEGR specially constructed with the quadratics of the torsion 
tensor [26]: 

T%, = - = e\ (dpet ~ duet). (2.3) 
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It is noticeable that the connection T* v = e\d u e^ used in TEGR not only describes 
a curvature-less geometry but also guarantees the independent parallel transformation 
among the four unit vectors e^x^), which form an orthonormal basis of the tangent space: 
• efl = r]AB-, where tjab = diag(l,—l,—l,—l) [25]. These vectors, e^, are commonly 
addressed by the components e A in the coordinate basis, i.e. = e A dp , while the metric 
tensor is given by the dual vierbein as g^ v = tjab £p . 

The field equations are obtained from the variation with respect to the vierbein, e4, 
which is the dynamical field of teleparallel gravity, given by 



fx 



\d^d^ - v{(f>) 



+4(V>(9A^) = e/, (2.4) 
where Qa u = e~ 1 5£ m /5e'^ denotes the energy-momentum tensor of the matter source, and 

G\ = 2e- 1 d,(ee A S/n - 2e\T% x S p ^ - \e\T 

is nothing but the equivalent geometrical structure of the Einstein tensor with Gp V = 
g^Gau, where 

with K^ u p = —^{r^ V p — T u ^ p — T p ^ v ^j . We can see that in the minimal limit (£ = 0), the 
field equation (2.4) is reduced to the simplest quintessence model. 

Torsion in (2.3) is defined by a Lorentz violated formalism for the tangent frame [26], 
which provides a technically simpler approach to acquire results of the cosmological interest 
[38]. The lack of Lorentz symmetry in teleparallelism results in extra degrees of freedom 
[33] as seen from (2.4), which in fact exhibits the equation of motion of all the 16 degrees of 
freedom of the vierbein field. To be more precise, given that the energy-momentum tensor 
Qp U and the Einstein tensor Gp U are both symmetric, the anti-symmetrization between 
indices A and v in (2.4) leads to a non-trivial constraint 

4^ (g^Sp^ - g^sM d x 4> = 0, (2.5) 



which forms a system of 6 equations in addition to the usual 10 equations from the sym- 
metric part of the field equations. Incidentally, this constraint vanishes if £ = so that 
TEGR shares the same number of dynamical degrees of freedom as GR. 

Nevertheless, we will see in the later discussion that only one extra scalar mode beyond 
the metric perturbation variables will be involved in the linear perturbations of teleparallel 
dark energy theories. However, this new degree of freedom is not dynamically independent 
due to the constraint in (2.5). 



1 We use the notations as follows: Greek indices n,v,... and capital Latin indices A, B,... run over 
all coordinate and tangent space-time 0, 1, 2, 3, while lower case Latin indices (from the middle of the 
alphabet) i, j, ... and lower case Latin indices (from the beginning of the alphabet) a,b,... run over spatial 
and tangent space coordinates 1, 2, 3, respectively. 
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The flat Friedmann-Robertson- Walker (FRW) background metric 



ds 2 = dt 2 - a 2 {t)5 i jdx l dx : > 



is conventionally resolved by the background vierbein choice 



4 



diag(l, a, a, a) 



(2.6) 



where a(t) is the scale factor. The Friedmann equations are given by 




) 



(2.7) 



(2.8) 



where 



P4> = -<A 2 + ^)-3^V, 

P<t> = \<t? ~ V(<f>) + 4£ff# + C {ZH 2 + 2H) <j) 2 



(2.10) 



(2.9) 



with H = a/a the Hubble parameter. The variation of the action with respect to the scalar 
field yields 



where Va, = dV/d(j). Consequently, the background equation motion of <f) is given by 



The continuity equation p + 3H(1 + w)p = holds for both p m and p^, where w = p/p 
is the equation of state parameter. It is interesting to note that, due to the lack of trivial 
Einstein frame available through conformal transformations, the effective equation of state 
of dark energy = p^/p^ in teleparallel cosmology can exhibit [29] quintessence-like, 
phantom-like or phantom-divide crosing evolutions, respectively. 

3 Perturbation Equations 

The general vierbein field has 16 components, which are dynamical in teleparallel gravity 
theories. In particular, 6 of the 16 degrees of freedom, arising as a result of the non- 
invariance of the local Lorentz symmetry, are unfamiliar in the metric perturbation sce- 
nario, and shall not appear in the field equation of the minimal TEGR limit. The perturbed 
vierbein field is given by [36] 



U<j> - £T(t> + = 0, 



(2.11) 



4> + 3H<P + QiH 2 ^ + = 0. 



(2.12) 



e° = <$(1 + + aS^diiF + a) + aS^G, + on) 
e« = a^(l -<p) + a5 i ^ 1 {d l d a B + d a C l + h a i) 
+ a5 i pL B a t + 5l{d a a + a a ) 1 



(3.1) 
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where B a i is an antisymmetric spatial tensor with B^j + Bji = and Bij = 5fB a j, while 
Ci, Gi, and ctj are transverse vectors with d l Ci = d l Gi = d l oti = 0. It is easy to obtain the 
corresponding metric from (3.1) as: 

#oo = 1 + 2tp 

g i0 = a(diF + G t ) 

gij = -a 2 [(l - 2ip)6ij + h i:j + d l d j B + djd + diCj], (3.2) 

which show the familiar metric perturbations in GR. Hence, we observe that the unfamiliar 
modes in the perturbed vierbein (3.1) are the scalar a, the vector a 1 and the tensor B a i, 
which comprise 1 + 2 + 3 = 6 degrees of freedom. Namely, the 16 components of the vierbein 
field are illustrated by the 5 scalar modes: tp, <fi, B, F and a, 3 vector modes: Q, Gi and 
«j, an antisymmetric tensor B a i, and a transverse traceless tensor mode hij, respectively. 

In the following study, we focus on the cosmological implication of the scalar per- 
turbations. To simplify the calculations, it is convenient to eliminate the variables F, B 
and Ci through the general coordinate transformations x^ — > x^ + e fl (x). The vierbein 
perturbation is provided by [36], given by 

e° = $(1 + V) + aS^dia 

= <(1 ~<P) + a^B a i + §d a a, (3.3) 

leading to the perturbed metric in terms of the longitudinal gauge form: 

ds 2 = (1 + 2ip)dt 2 - a 2 (l - 2ip)5ijdx l dx j . 

Subsquently, we obtain the torsion tensor and scalar to be 

T° i = -ditp + ad dia, 

Toj = (H - 0)<5* + cWj - a^dj&a, 

T l jk = d k {8)v - B*j) - djiSicp - B\), (3.4) 

and 

T = -6H 2 + l2H(tp + HiP) + Aa- 2 d 2 a m , (3.5) 

respectively, where a m = aHa. 

On the other hand, we can also decompose the scalar field to homogeneous and per- 
turbed parts: <j) — > <p(t) + 8(j)(t,x l ). The scalar perturbations of the matter source Q„ are 
denoted by the perfect fluid perturbations: 0q = — (p m + Sp), Gf = — (p m + p m )di5u and 
©j = (P m + fip)fiji wnere $ u is t ne scalar vector potential of the fluid. 

After substituting these expressions to the field equation (2.4), we derive the pertur- 
bation equations given by the energy density ©o°: 

k 2 



(1 + k'&* 



6H(<p + Hi/>)-2-~<p 





a z 



+ k 2 (y+SQ - <j> 2 ^J - 6CH 2 (f>5^j = -k 2 5 P , (3.6) 
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the energy flux ©o 4 : 

2(1 + k 2 ^ 2 )^ + Hip) - K 2 4>d l b4> + Ak 2 ^ \&ip + 
= -n 2 ( Pm +p m )d i 5u, (3.7) 
and the momentum density ©.;°: 

- 2(1 + K 2 £4?)di{ip + Hip) + K 2 <j>di5<p + 4K 2 ZH<pdi5<P = K 2 {p m + p m )di6u, (3.8) 
respectively, while the momentum flux ©^ is divided by its diagonal components (i = j): 
(1 + k 2 ^ 2 ) \l2H((p + Hip) + 2(<p + Hip + 2Hip) 

-k 2 - - 2 j + 4k 2 £<^>(<^ + - ±K 2 £H4>5(p = K 2 5p , (3.9) 

as well as its anisotropic parts (i / j): 

2k 2 
£T(1 + k 



Similarly, the perturbed equation of motion for the scalar field (2.11) is given by 

5<p + 3H50 + + ^ + 6£# 2 (^ - 2(<£ + 3il<^ - ^(30 - ^) 

A: 2 

-l2^(ip + HiP)-^cP^a m = 0. (3.11) 

To complete the perturbation equations of teleparallel dark energy models, we shall 
also include equations from the constraint (2.5). In linear perturbations, this constraint 
yields 

00V + ^4>djB ji + Hd i S(p = 0, (3.12) 

which is necessary for the consistent relation of the matter source ©o« = ©io, as also 
indicated from the comparison between (3.7) and (3.8). Given that the antisymmetric 
tensor Bij satisfies didjB 13 = 0, we find that the gradient of (3.12) leads to 

ip = — ^5<p. (3.13) 

<P 

Note that this equation holds even for the absence of the matter source with the 
non-minimal coupling £ 7^ of the theory (2.1). 

4 Matter Growth 

Although teleparallel gravity gives a gravitational concept different from GR, the perturba- 
tions of the matter sector shall be treated in the usual manner regardless to the geometrical 
background. In particular, the equation motion of a particle in the presence of gravity in 
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TEGR is mathematically equivalent to the geodesic equation in GR [26], while the conser- 
vation of energy-momentum V^0^ = is the covariant derivative with respect to the GR 
connection [33, 39]. Thus, in teleparallel cosmology, the gauge-invariant matter density 
perturbation of the dust-like source, 5 m = 5p m /p m + 3H5u, also satisfies the perturbation 
equation 

S m + 2HS m - 4TTG eS p m S m = 0, (4.1) 

where the effective gravitational coupling constant G e s coincides to the Newton's constant 
in the minimal case of £ = 0. In the discussion, we consider the Fourier space by simply 
replacing the operator d 2 with the wave number k 2 . In the sub- horizon regime k 3> aH, 
the perturbed equation of V^G^ = Cctn lead, to h. useful relation 47r(j e ff Pm$m 

~ k 2 ip/a 2 

under the longitudinal gauge. The effective gravitational coupling G e g can be obtained 
once the relation between ip and 8 m is found. 

To proceed a further simplification, we adopt the quasi-static approximation 

\X\<\HX\, for X = (f>,ip,tp, 

to the perturbation equations (3.6)-(3.11). Together with the sub-horizon condition k S> 
aH, we find that (3.11) becomes 

% + V H ) 84> - 4£<^a m ~ 0, (4.2) 



while (3.6) gives 



k 2 

2(1 + K V)-V + ^W--^. ( 4 - 3 ) 



a 2 

where |£| is O(10 -1 ) according to the observational constraints [32]. We shall further 
neglect the contributions of V$ and as for the general scalar-tensor theory of dark 
energy [40], and it can be checked that any potential for teleparallel dark energy indeed 
satisfies k 2 3> a 2 KV$ and k 2 3> a 2 V^ in the sub-horizon regime [32]. 

From (3.13) and (4.2), we obtain that 5<j) = -<j>ip/H ~ 4£0a m . The relation (3.10) 
results in 

*-{? + t&&)* (44) 

where e = K 2 <p 2 /2H 2 . From (4.3) and (4.4), we finally arrive at 



which implies 

G_ 

1 + K 2 £4P ) 1 +~K?£(j) 2 

Note that e is proportional to the kinetic energy density of the scalar field, which is sub- 
dominated to the energy density of dark energy at the present in teleparallel dark 
energy models. In the minimal quintessence limit £ = 0, the cosmic acceleration in fact 
requires e <C and the effective gravitational coupling reduces to the expected result, i.e. 
G cff = (1 + e)G ~ G. 
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4.1 Purely gravity coupling V = 

In the following, we exam the matter growth in teleparallel dark energy models with specific 
forms of potentials. The simplest tracker behavior has been found in the potential-less 
case (V = 0) where one can achieve the late-time cosmic acceleration for £ < 0, while 
the observational constraints from type la supernova (SNIa), baryon acoustic oscillation 
(BAO) and cosmic microwave background (CMB) data reveal a theoretically viable region 
of — 1 < £ < —0.125 with the best-fit value at £ = —0.35 [37]. To analyze the evolution in 
this case, we re-parametrize two new variables: x = K<j)/y/6H and y = \/—£,K<p so that the 
Friedmann equation (2.7) is rewritten as 

l = n m + x 2 + y 2 . (4.7) 

The evolutions of x and y with respect to N = In a are given by 

dy 



iN V^|- (4-8) 
Note that (2.8) also provides a useful relation: 

(l-y 2 )§2=-3x 2 + 2v^xy-^n m , (4.9) 

to solve the evolutions numerically from some given initial conditions. 

By using Q = d\a.{5 / a) / d\a. a to re-parametrize the matter perturbations, the second 
derivative equation (4.1) is reduced to 

w + ( 2 + w) {G + 1) + {G + 1)2 = \ Q n ™ ' (410) 

where Q = G e f[/G. Subsequently, from (4.6) we find that 

We remark that (4.10) is applicable to the generic modified gravity theories if the effec- 
tive coupling G c ff is given. In the GR case, the contribution of the Q 2 term in (4.10) is 
negligible throughout the growth history from matter domination [41]. Here, we adopt 
the same approximation given that the present value of Q is not smaller than —0.5 for all 
of the considered cases in this work. As a result, we can have the evolution of Q solved 
simultaneously with the other parameters x and y from given initial conditions. In Fig. 1, 
we show the evolution of the growth rate 5 as a function of the redshift z in the purely 
gravity coupling teleparallel dark energy model (V = 0), where the dotted, red-solid and 
dot-dashed lines correspond to the observational lower bound, best-fit and upper bound 
values of £ = —1, —0.35 and —0.125 respectively, while the green-solid line represents 



-8- 



0.2 







0.90 
















0.85 










0.80 










loO.75 










0.70 










0.65 








0.60 










v, 0.0 0.1 0.2 0.3 0.4 0.5 








Z 



1 2 3 4 5 6 

Z 



Figure 1. Growth rate S as a function of the redshift z in the purely gravity coupling teleparallel 
dark energy model (V = 0), where the dotted, red-solid and dot-dashed lines correspond to the 
lower, best-fit and upper values of £ = — 1, —0.35 and —0.125, respectively, while the green-solid 
line represents £ = —0.35 with a fixed value of Q = 1. 



£ = —0.35 with a fixed value of Q = 1 so that gravitational interaction mimics that of 
GR. The initial values provided are {tt m = 0.999 , y = 1.0 X 10~ 6 , <7 = 0} far z » 1. The 
results in Fig. 1 indicate that the gravitational interaction is in fact enhanced (Q > 1) 
during the evolution, which ensures a more reasonable approximation by dropping the Q 2 
term in (4.10). 

4.2 Exponential potential V = Vq e~ 2KX ^ 

We can extend the scenario to teleparallel dark energy models with non-vanished potentials. 
As an illustration, we concentrate on the exponential type, V = Voe _2KA< ^, which has a 
constant potential limit V = Vo if A = 0. We shall introduce a new parameter, v = 
K\fV /\/3H, to address the evolution of the potential. The Friedmann equation (2.7) is 
rewritten as 

l = n m + x 2 + y 2 + v 2 , (4.12) 

where x and y share the same definitions as the previous case. The evolution equations of 
these components are given by 

dx / H 



d y r~Fi 
In = V ~ 6 ^ 

+ (4.13) 

We consider the regime of £ < for the exponential potential since it is favored from 
SNIa, BAO and CMB observations [32] and solve x, y and v together with the density 
perturbation Q = d\n(5 / a) / d\n a in (4.10). Note that (4.9) is unchanged. 

In the case with the non-minimal coupling of £ < —0.01, the result (refer to the green- 
solid line in Fig. 2) is similar to the quintessence model where the potential dominates 
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Figure 2. Growth rate 8 as a function of z in the teleparallel dark energy model with V — Vq e~ 2/tA ^, 
where the red-solid, blue-solid and green-solid lines correspond to £ = —0.4, —1 and —0.125 with 
the initial values of {£l m = 0.999 , v = 1.0 x 10~ 6 }, and the dotted line depicts £ = —0.4 with the 
initial values {fi m = 0.990 , v — 6.0 x 10~ 3 }, while we have fixed A = 0.5 and the initial values for 
z > 1 to be {y = 1.0 x 10~ 6 , Q = 0}. 



during the evolution. However, if £ is around the best-fit value of ~ —0.4 [32], the purely 
gravity coupling tracker behavior appears when the given initial value of v is sufficiently 
small, as illustrated by the dotted line in Fig. 2 with the initial value v ~ irr 6 . The 
comparison to the red-solid line shows that the gravitational interaction is weaker if the 
cosmic acceleration is driven by the potential instead of the NGF tracker effect. On the 
other hand, if the non- minimal coupling £ is ~ —1 as described by the blue-solid line in 
Fig. 2, then only NGF tracker behavior is found regardless to the initial values of the 
potential. Clearly, the parameter y will finally dominate the dark energy density $7^ at the 
present. 

5 Conclusions 

In this work, we have applied the vierbein perturbation scenario to study the cosmological 
perturbations in teleparallel dark energy models with the non-minimal coupling between 
the scalar field and gravity. Unfamiliar with the metric perturbations, a new scalar mode, 
a m , appears in the linear equations and mainly contributes to the anisotropic between the 
gravitational potentials ip and (p. This new scalar degree of freedom is fully suppressed in 
the sub-horizon scales of f(T) gravity [36], but has significant contributions in the large 
scale structures [35]. Nevertheless, the dynamical degrees of freedom of the theory do not 
increased by a m given that the constraint in (2.5) provides an additional algebraic relation 
to the scalar mode variables. 

To study the matter growth in teleparallel cosmology, we have used the quasi-static 
approximation in the sub-horizon scales. The effective gravitational coupling becomes scale 
independent in the sub-horizon regime similar to the generic scalar-tensor theories. In the 
case where the cosmic acceleration is driven by the non-minimal gravity-field coupling term, 
i.e. the £i7 2 </> 2 term in (2.9), we have found that the gravitational interaction is enhanced 
so that matter grows faster than that in GR. In the broader case where the scalar field 
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has a potential, the quintessence behavior recovers with a small coupling limit of £, while 
the non-minimal gravity- field tracker is followed when £ ~ — 1. In the case with £ around 
the best-fit value (~ —0.4) from observations, we can also observe a faster growth rate if 
the NGF tracker is followed instead of the usual quintessence-like evolution. These results 
feature the matter perturbations in teleparallel dark energy models and worth for further 
investigations to look for the possible deviations from GR. 
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